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ABSTRACT 


The sea breeze circulation was investigated using a combination of acoustic doppler 
sodar, doppler lidar and conventional observations in central California on the coast of 
Monterev Bay in September 1987. The study was called LASBEX (Land Sea Breeze 
EXperiment) and used the combined effort and resources of the Naval Postgraduate 
School (NPS), NOAA Wave Propagation Lab (WPL) and Naval Environmental Pre- 
miction Research Facility (NEPRF). 

A monostatic three-axis phased-array doppler sodar was able to detect the various 
features of the sca breeze front, which has many of the characteristics of a gravity head 
current produced in laboratory experiments, such as the longitudinal vortex, lobe and 
cleft structure, and strong regions of turbulence and convergence. Profiles of the wind 
field and its component variances Were generated every five minutes. Vertical velocities 
ШРОТ О ms were detected in the updraft region at the front. Vertical velocity fields 
were used to derive the horizontal divergence associated with the sea breeze and values 
Or convergence in the range of J x 10-*s-'! and greater were calculated, These are believed 
to be more accurate values than anv found in previous sea breeze studies. 

Doppler lidar radial velocity information was useful in mapping the various lavers 
плашта structure of the sea breeze. The wind field of the sea breeze seen by the lidar 
was compared with the sodar. The results of this comparison showed that the sodar was 
well suited to investigate features under 600 meters in height, such as the updraft region 
behind the front and the variance of the wind, whereas the lidar was better at describing” 
the larger scale flow, such as the return flow of the sea breeze. the convergence zone 
along the sea breeze front and prevailing svnoptic winds. Together, these instruments 
gave a comprehensive and complete description of the sea breeze frontal structure, evo- 


lution and movement. 
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І. INTRODUCTION 


A. GENERAL DESCRIPTION 

The land-sca breeze is a diurnal mesoscale circulation that occurs on sea coasts due 
to the differences in surface temperatures between the land and sea. À typical horizontal 
temperature gradient is 1° C;20 kin (Atkinson, 1981). This situation induces a thermally 
direct circulation that alternates between onshore or sea breeze and offshore or land 
breeze. The circulation occurs during the day due to solar heating of the land which 
warms faster that the sea. This results in a horizontal pressure gradient with slightly 
higher pressure over the land. The pressure differential causes air to flow at higher levels 
towards the sea, thereby leading to convergence aloft. The resulting counter-flow at the 
surface is the sea breeze. Fig. 1. shows a lake breeze circulation, which has a similar 
structure as the sea breeze circulation. 

ше заа breeze is due to the sane thermal structure except it is reversed. At night, 
the land loses heat quicker than the sea and, therefore, sets up a circulation opposite to 
Eiessea breeze. In both cases, the flow at the surface is from the cooler region toward 
the warmer region. By comparison, the sea breeze is the stronger circulation. since the 
heating of the land provides for a more intense heat source to drive the flow. 

The sea breeze occurs more frequently given a weak horizontal pressure field 
meetojue, 1962: Simpson ct al., 1977). This is because the thermal eflects have to alter 
anv eXisting conditions. If the existing pressure field is strong enough, the resulting 
gradient winds can overwelm the land-sea breeze effect. 

Under some conditions, an offshore gradient wind will generate a frontal zone at 
Ме сасе сасе of a sea breeze (Atkinson, 1981). Estoque (1962) notes in his model 
that un offshore flow counter to the sea breeze will enhance the horizontal temperature 
gradient. The sea breeze front is also characterized by gradients in humidity, wind speed 
and wind direction. This sea breeze front is similar in form to a synoptic scale cold front, 
because the cooler moist marine air is advancing into the warmer drv land regime. 
Passage of such a front is a measurable feature and can induce significant mesoscale 
meteorological phenomena such as showers, thunderstorms and gusty winds (Pielke, 
1973; John, 1980; Richiardone and Pearson, 1983). The horizontal extent of the sea 
breeze is large enough that it can develop a component parallel to the coast due to 
coriolis deflection (Defant, 1951; Walsh, 1974). 
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Но. 1. Typical Lake Breeze Circulation. (from Кееп апа Lvons, 1978). 


B. MOTIVATION FOR THE STUDY OF LAND-SEA BREEZE 

The land-sea breeze circulation and sea breeze front can affect the dispersion of 
aerosols and pollutants over a wide area. This effect has been observed and modelled 
many times (Lyons, 1972; Lyons and Olsson, 1973, Anthes, 1978; Kondo and Gay 
1979; Kitada et al., 1986: Ogawa et al., 1986). In a sea breeze, particulates and gases 
can be carried inland over long distances. Collis and Russell (1976) summarized some 
of the important effects that particle pollutants have on climate, such as cloud forma- 
поп, precipitation and transfer of solar radition. Pollutants can enhance cloud forma- 
tion and result in acidic rainfall. Particulates and gases can block incoming solar 
radiation and trap outgoing long-wave radiation. 

In the case of lidar (light detection and ranging) sensors, concentrations of 
pollutants such as .VO, and SO, can attenuate or distort the return signal so that infor- 


mation may be misinterpreted (Measures, 1984). 


Measuring the intensity and mapping the structure of the sea breeze can give a 
better understanding of how particles are transported and dispersed. Satellite data could 
be used to detect a sea breeze occurrence. In addition, forecasts of pollution events may 
be feasible. Surface data also can be used to evaluate satellite observations of the 


planetary boundary layer and the mesoscale features that occur during the sea breeze. 


C. PREVIOUS STUDIES 

A study was conducted by Johnson and O’Brien (1973) of the sea breeze along the 
Oregon coast. In their study, the sea breeze was found to have a maximum depth of 700 
ВЕЕТ, With a maximum speed of 7.5 m's at the surface. The zonal wind was utilized 
to describe the variability of the sea breeze because of the north-south orientation of the 
Oregon coast. The imlai.d penetration of the sea breeze front was up to 60 km. A return 
flow was found between 1060 and 1500 meters. Johnson and O'Brien made the following 
observations: the low-level onshore flow of the land-sea circulation was located entirely 
Metin the marine layer; the marine layer tends to deepen briefly at the onset of sea 
eae due to the low level convergence ahead of the sea breeze front; the return flow 
above the inversion appears to be in the form of surges, probably in response to the 
surges observed in the sea breeze; and a distinct wind maximum follows the sea breeze 
front inland. 

Simpson, et al. (1977) observed sea-breeze fronts that penetrated over 45 km ınland 
with some fronts moving as far as 80 km. These fronts had an average advance speed 
Eg ons They also showed that the structure of the sea breeze front is similar m 
structure to gravity currents produced in laboratory water-tanks. 

Mitsumoto et al. (1983) conducted a laboratory experiment which used a water-tank 
to simulate the sea breeze. This experiment reproduced many of the characteristics of 
the sea breeze flow, such as strong vertical motion at the leading edge of the sea breeze. 
compensating flow aloft and production of longitudinal convective vortices within the 
a breeze flow., 

Ogawa ct al. (1986) observed the gravity head structure reproduced by Mitsumoto’s 
laboratory experiment during a lake breeze. Ogawa's experiment utlized a kytoon- 
mounted ultrasonic anemometer-thermometer to measure instantaneous wind and tem- 
perature in order to observe the turbulent structure of the lake breeze. He found that 
the penetration of a lake breeze produced an upward rolling motion within the head of 
the front. He recorded upward flow velocities around 3 m's. Ogawa found that the 


temperature profile had peaks indicating the frontal passage. These peaks were due to 


the boundary between the offshore wind aloft and the lake air below. Another feature 
that was observed was the development of the Thermal Internal Boundary Laver 
(TIBL). The TIBL develops due to the temperature discontinuity between the land and 
the water (Stunder and Sethurman, 1985). The TIBL plays a significant role in the 
fumigation of coastal pollution that 1s advected into it. 

Estoque's (1962) model had a maximum penetration of sea-breeze front of 32 km 
and a minimum penetration of 18 km. The maximum penetration occurred during calm 
initial conditions and the minimum penetration occurred during offshore counter-flow. 

Pearson et al. (1983) used a model and found that with a mean flow perpendicular 
to the coast, the speed of the sea breeze front is a linear function of the onshore com- 
ponent of the mean flow. He found that the maximum vertical velocity in the front at 
any given time was not affected bv the mean flow. Another characteristic that Pearson 
found was that in a negative shear case (upper flow is from land to sea), there was an 
increase im the vertical extent of the sea-breeze front. 

For a general descriptive summary of sea breeze studies, Atkinson (1981) provides 


a comprehensive review from many locations and experiments. 


D. LAND-SEA BREI ZE EXPERIMENT (LASBEX) 

A measurement program was conducted in the Monterey Bav area in September 
1987, to understand the structure of the land-sea breeze effect and generate data to 
evaluate many of the models which have been used to demonstrate the behavior of the 
air flow. This experiment was unique because it utilized a varied array of instruments 
including doppler acoustic sounder (sodar), a doppler lidar, rawinsonde launches over 
the the sea, radiosonde launches over land, conventional surface recordings of temper- 
ature, hunudity and wind and corresponding observations by the NOAA 9,10 and GOLS 
satellites, This combined instrument analysis allowed for comparison of the data among 
many sources. The final result will be a comprehensive view of the land-sea breeze in- 
tensity and structure. This study focuses on data from one sodar supplemented by hdar 
observations. 

The products obtained from the sodar data were the height of the marine boundary 
layer and time-height cross sections of wind direction, wind speed, vertical velocity and 
the standard deviation of the vertical velocity. From these it was possible to compute 
time-height cross sections of horizontal divergence and turbulent kinetic energy. The 


A 


tine-height cross sections will be compared with conventional data obtained from 


radiosonde, rawinsonde and surface observations. The boundary laver height observa- 
tions will be compared with the time-height and conventional data. 

The lidar provided pictures of the radial wind field in the sea breeze front. The lidar 
observations will be compared with the sodar to find differences between the instruments 
ability to map the sea breeze front and measure the associated wind field. 

The combination of lidar and sodar has been applied before to analyzing the at- 
mosphere. However, the two instruments have not been used to study the sea breeze 
circulation. Russell et al. (1974) used the combination to observe the atmospheric 
boundary layer in an urban environment. One comparison in Russell’s study was the 
ability of either instrument to detect subsidence. In the case of the lidar, detection was 
degraded because aerosols tend to be excluded from descending air. For the sodar, the 
temperature gradient associated with a large-scale subsidence inversion is easily detect- 
able. Lopez (1977) used the combination to study the properties of convective activity 
and small cumulus formation. This combination is effective since either instrument can 
detect different features to give a better overall picture of the atmosphere. 

Y The ОБЕ О: IMIS thesis 15 to evaluate the eflectiveness of the sodar in the de- 
tection and mapping of sea breeze events and to compare its measured wind field with 
that of the lidar’ A case study on 16 September 1987 is presented. An overall look at 
tvpical characteristics cf the sea breeze in the Monterey Bay. Salinas Valley will be gen- 
erated for use in other studies, such as aerosol transport or the evaluation of satellite 


imagery. 


П. DESCRIPTION OF EXPERIMENT 


A. LOCATION OF THE EXPERIMENT | 

LASBEN was conducted from between 15 and 30 September 1987 on the Monterey 
Вау coast. The collection of data started at 1200 Pacific Standard Time (PST) on 15 
September and ended at 1300 PST on 30 September. All instruments were operated 
continuously except the contracted research vessel Silver Prince, which was stationed in 
Monterey Bay to monitor surface wind, temperature, humidity and pressure as well as 
launch rawinsondes. The ship collected data from 1000 PST until 1600 PST on various 
davs throughout the experiment. The instruments evaluated in this paper were the 
Pennvlvania State University (PSU) sodar and the NOAA Wave Propagtion Laboratory 
(WPL) lidar. The sodar was located southeast of Castroville at the mouth of the Salinas 
Valles. The lidar was located on a flat knoll approximately 1.5 km from the coast, 
southeast of Moss Landing. It was on a line of 310° with the sodar, which approximated 
the orientation of the Salinas Valley. The separation between the sodar and the lidar 
was 5.5 km. The Silver Prince was stationed in Monterey Bay approximately 16 km from 
the lidar location. The three main instrument locations are shown in Fig. 2. 

The terrain around the sodar site was at elevations of less than 30 meters, with ele- 
vations greater than 30 meters northeast, starting 2.5 km away. The interference caused 
by the town of Castroville was considered minimal for the experiment. The terrain sur- 
rounding the lidar was generally flat with higher elevations to the north and east. To the 
east the terrain rose to levels of 100 meters and greater, up to 400 meters. To the north, 
the terrain generally had a maximum elevation of 150 meters. North-northeast of the 
lidar, Elkhorn Slough, a brackish inlet, interruped the elevated terrain. 

At the lidar, sodar and shin, twenty-second averages of vector wind velocity, air 
temperature, relative humidity and surface pressure were collected at the stations during 
the experiment. Radiosondes were launched at the sodar station at the discretion of the 
research team, coinciding with significant changes in the environment. 

Satellite AVHRR imagery during LASBEX was available from the Scripps Institu- 


tion of Oceanography. 


B. ACOUSTIC DOPPLER SODAR 
During LASBEX, the doppler sodar was used to monitor the acoustic backscattered 


energy of turbulent fluctuations in the atmosphere. The sodar is a 3-axis monostatic 
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Fig. 2. Location of Silver Prince, Lidar and Sodar: Тһе Silver Prince was at lo- 
cation A. The Lidar was at location B. The Sodar was at location C. Ele- 


vations are in feet. Distances are in kilometers. 


doppler system inanufactured by Xontech and consists of a phased array of 25 
transducers which operated at a frequency of 1600 hertz. The pulse repetition interval 
is 10 seconds per beam, giving a cycle time of 30 seconds between the 3 acoustic beams. 
Two cvcles were integrated per nunute. A wind profile was sampled every five minutes, 
giving an effective number of 10 cycles per record. 

The sodar uses sound energy to detect the small-scale temperature fluctuations and 
discontinuities in the mean temperature profile (Crease et.al., 1977; Koracin and 
Berkowicz, 1988). The addition of doppler capabilities allows for the recording of wind 
speed and direction (Kaimal and Haugen, 1977), The backscattered energy provides in- 
formation about the height of the marine boundary laver, wind speed, wind direction, 
vertical velocity and vector winds. Neff and King (1988) demonstrated that the sodar 
can give indirect spatial and temporal information on both the location and thickness 
of turbulent layers. In addition, Neff was able to identify nuxing processes in the 
boundary laver. Small mixing processes could be attributed to turbulence, where as 


large processes could be associated with effects such as Kevin-Helmholtz instabilities. 


C. DOPPLER LIDAR 

The lidar used during LASBEX was the WPL 10.6 micron doppler systemy a mmama 
operated under computer control and had real-time data processing and color display 
capabilites. The system could compute the zero, first and second moment of the echo 
spectrum for each range gate. For LASBEX, the lidar was operated in the range-height 
indicator (RHI) and plan postion indicator (PPI) modes. In RHI, the beam of the hdar 
sweeps in a 180? arc from foward to overhead and then behind the instrument. In PPI. 
the beam sweeps at a fixed angle of elevation from the horizon in a 360* circle around 
the lidar, the path describing an inverted cone. In both modes. the display of informa- 
tion is of the radial velocity relative to the lidar. Famuharity with a doppler presentation 
Is required to analyze the PPI display. This analysis 1s similar to a doppler radar PPI 
and is described in Wood and Brown (1983). In RHI, the hdar beam was aligned parallel 
to the axis of the Salinas Valley with an orientation angle of 140°, or east-west with an 
angle of 90°, with the preferred orientation being east-west. This was because 
climatological records indicated that the sea breeze blew from due west in the Moss 
Landing area. In PPI, the elevation of the lidar beam was 1° and 5°. The 1° beam 
shows the low level structure below 200 meters, whereas the 5° beam gives structure 
above 300 meters. These values were determined triginometrically. Table 1 lists some 


of the clear-air radial velocity parameters of this instrument: 


Table 1. WPL DOPPLER RADAR CAPABILITIES 


Accuracy of radial velocity measurment 


Operation modes PPI, VAD, RHI, range-time. raster-scan 













The lidar uses the backscattered energy of a laser from gas molecules and 
particulates to detect atmospheric features such as clouds or aerosol layers. Collis and 
Russell (1976) cited an example of how lidar was used to find the depth of the mixed 
laver or determine the heights of inversions, which can affect the dispersion of 
pollutants. Noonkester ct al. (1974) was able to detect Kelvin-Helmholtz instabilities in 
the boundary-layer. In clear air, a lidar can detect the stratification associated with a 


deep stable boundary layer. 


D. SURFACE MEAUREMENTS 

A portable coastal meteorological station was used to take 20-second averages of 
wind direction, wind speed, air temperature, relative humidity and barometric pressure 
at the surface. This instrument used a propeller anemometer to measure wind speed. 
A self-aspirated thermometer and hygrometer measured air temperature and relative 
humidity. Wind direction was taken from a vane that gave angles with respect to true 


ШІ. Lhe ranges of the instrument errors were: + 5° for wind direction: + 0.3 m s for 


wind speed; + 0.3° for air temperature: and + 5%o for relative humidity. 


E. SYNOPTIC SITUATION BETWEEN 16-30 SEPTEMBER 1987 

On 16 September, California was under the influence of a surface high pressure 
svstem olf the coast aud a local thermal trough in the center of the state. The winds 
were northerly to northwesterly between 5 and 10 knots. This configuration continued 
to influence the west coast from 17 to 24 September. Fig. 3shows the synoptic situation 
during this period. Around the 25th, a short wave which was detectable at the $50 mb 
level began to pass through northern California from the northwest. This wave built up 
a strong thermal gradient off the coast of California, providing an increase of cold ther- 


mal advection. This was the most significant feature associated with an actual synoptic 


frontal passage during LASBEX. Fig. 4 shows the short wave location. Winds during 
this Wave passage were mainly northerly with speeds averaging 7 m's or more. After the 
passage of the wave on the 26th, high pressure off the coast and the thermal trough in 
the center of the state were reestablished. This weather pattern persisted for the re- 
mainder of the experiment. 

The combination o a high pressure system over the eastern North Pacific Ocean 
and a thermal trough in central California is the typical summmer and early autumn 
synoptic pattern. This situation accounts for the common Monterey Bay pattern of low 
clouds and fog with light and variable winds during the night and morning. A tropical 
disturbance named ‘Otis’ formed west of the Baja Peninsula during the 23rd but did not 


significantly affect surface and lower level flow during that part of the experiment. 
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Fig. 3. Typical Surface Synoptic Situation 17-30 September 1987: Solid lines are 


sea level ısobars. Values are + 1000 mb. ~ 


25 September 





Fig. 4. 850mb Synoptic Situation 25 September 1987: Solid lines are isoheights. 
Values are + 1000 meters. Dashed lines are isotherms. Temperatures аге 


in degrees Celsius. 


Ш. DATA AND ANALYSIS 


A. SODAR THEORY 


l. Basic Sodar Theory and Turbulence in the Inertial Subrange 


This discussion of sodar theory and operation follows Underwood (1981). The 


sodar utilizes the backscatter of acoustic energy to determine certain parameters in the 


atmosphere. The relatronship between the acoustic power emitted and the backscatter- 


ing cross section is defined by the sodar equation. 


7 = | 4, 
P(R)- P, p Le Кы к ye 


an БЕ Ar 
where 
P(R) 2 received acoustic power (watts), 
P, — transmitted acoustic power. (watts), 
4 = cross sectional area of the receiver antenna (m), 
G — antenna gain factor. 


L, — length of the acoustic pulse in space (m). 


ЕЕ. , : = 
« = path averaged molecular attenuation coefficient (m ), 


КЕ Zdıstanee from the scatterer to the receiver (1m). 
= “excess” attenuation and 


o(z) = backscattering cross section (m™'). 


The backscattering cross section can be related to the temperature structure 


parameter C?, and the velocity structure parameter C?, as shown by Little (1969). These 


structure parameters indicate the intensity of high frequency variations of temperature 


or velocity. Explicitiy, the backscatter cross-section is 
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wliere 

o(z) = backscattered acoustic energy, 

с = speed of sound, 

T = absolute air temperature, 
2 - 

C, — velocity structure parameter and 
2 

C7 = temperature structure parameter. 


Turbulent mixing m clear air causes fluctuations in potential temperature and 
Water Vapor Which scatter and attenuate the acoustic energy in the sodar signal. The 
backscattering cross section determines how much energy is returned to the sodar an- 
tenna after attenuation by the turbulence. Turbulent eddies with a size equal to one half 
the wavelength of the emitted acoustic energy Aas (Russell et al., 1974) ате аепев је 
with a sodar. Eddies of this size fall within the inertial subrange of turbulence and are 
assumed to be isotropic. Energy relationships in the inertial subrange follow the 
Kolmogoroff description of isotropic turbulence. The backscattered data are processed 
using Fast Fourier Transforms (FFT) to establish the frequency shift due to the radial 
components. The two main fields which the monostatic doppler sodar can map are the 
thermal structure and the niean velocity field along with its variances. 
а. Thermal Structure 
Molecular attenuation is described by Sutherland (1975) and Piercy et al. 
(1977) as those effects at the molecular scale which influence the propagation of energy 
through a medium such as air. The factors which cause molecular attenuation of sound 
energy are due to thermal, viscous and diffusive effects as well as rotational and 
vibrational losses on the molecular scale. The molecular attenuation increases with fre- 
quency and ts the limiting factor in the maximum usable frequency in a sodar. The lower 
limit is based on antenr 1 size. “Excess” attenuation is due to turbulent spreading of the 
beam pattern (Brown and Clifford, 1976). Any of the sodar signal that is not attenuated 


can be scattered back to the receiver by turbulent eddies according to (3.2). 
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As the lapse rate changes in the atmosphere, any mechanical turbulence will 
create small eddies that backscatter the acoustic energy. The monostatic sodar (0 = 
150°) gives an indication of the changes in the lapse rate along the beam path. which can 
em lated to thermal structure (WyckolT et al., 1973). In a monostatic sodar, only 
temperature fluctuations contribute to received power according to (3.2). 

b. Velocity Structure 

Doppler processing is used to derive the velocity structure. This discussion 
follows Peterson (1988) in a handbook on doppler wind profilers. 

(1) Basic Setup of A Doppler Sodar. The sodar emits sound energy at a 
narrow band of frequencies and is reflected by turbulent eddies and temperature fluctu- 
ations. A frequency shift (doppler shift) 1s caused by any motion on the same axis as 
the sodar beam and is proportional to the speed of that motion. By measuring the 
doppler shift, the radial velocity of the scattering volume can be calculated. In order to 
obtain a wind velocity, Doppler shifts from three axes are required. The configuration 
of the sodar is to have one beam directed overhead, with the other two at perpendicular 
cardinal headings (ic. north and east). The beams are typically pointed at high ele- 
vations, so that they are sampling nearly the same flow. The vertical velocity from the 
beam directed upwards can be subtracted from the other two, leaving the horizontal 
components of the wind vector. 

The sodar used during LASBEX consists of 25 separate transducers. 
Using electronic phasing, these transducers Were fired at different times in order to direct 
the acoustic beam. If there was no phasing between all the transducers. the beam would 
travel directly overhead. Bv delaving the time that each transducer emitted. the beam 
of acoustic energy can be directed at different angles. 

(2, Range Aliasing. The signal from the sodar must be of a pulse гере- 
tition frequency (PRF) that precludes aliasing. Aliasing occurs when the the signal from 
the current pulse is overlapped by return signals from the previous pulse. These over- 
lapping signals will be erroneously interpreted. This tvpe of aliasing is called range ali- 
asing becuuse a return from a distant scatterer could be mistaken for a closer one. The 
time between pulses must be long enough to prevent this effect. The effective height that 


the sodar can function without range aliasing is determined by the following equation: 
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where 


Rmax = Maximum range, (n) 
c — speed of sound in air (m/s)and 

T =] 
PRF= pulse repetition frequency. (ame) 


The maximum height that the PSU sodar was set to detect is 750 meters, which is about 
half its maximum range, and only information starting 50 meters above the surface was 
recorded. These settings effectively avoided anv range aliasing during the experiment. 

(3, Velocity Aliasing. Only doppler shifts of certain frequencies can be 
determined without error. Doppler shifts that are either too large or too small will be 
incorrectly mistaken for other velocities. This is called velocity aliasing. The frequency 
limits for which a doppler shift can be interpreted correctly rely upon the sampling rate. 
The optimum selection is to have a sampling rate as high as possible. However, there 
must be balance between velocity aliasing which requires a high PRF and range aliasing 
which requires a small PRE. The wind speed data during LASBEX showed no extreme 
values and therefore velocity aliasing was negligible. 

/4, Data Processing. For each detection cycle, many samples are gath- 
ered and converted into a frequency spectrum using the Fast Fourier Transform tech- 
nique. Doppler shifts in a frequency spectrum show up as peaks offset from O or the 
no shift position. Motion towards the sodar 1s indicated by a positive doppler shift. A 
negative shift indicates motion away from the receiver. The doppler shift 1s converted 
into a radial velocity component of the wind. The height of the frequency peak is a 
measure of received power and the width of tlie peak a measure of spread of the acoustic 
energy over frequency. The spread of energy is directly related to the standard deviation 
of radial velocity. Using deviations of radial velocities among the three beams. the 
standard deviations of each velocity component (o ¢,, ¢,) can be calculated. 

The preferred method is to sample in the time domain, so that several 
consecutive samples can be gathered and averaged before conversion into the frequency 
domain with an FFT. A mean frequency spectrum is calculated from several consecutive 
frequency domains to generate a spectral average, which gives a better indication of a 
doppler shift in the sigual. The background noise is filtered out to better refine the fre- 
quency peak. This cycle of sampling, conversion, averaging and filtering 1s set bv the 
user and typically lasts one or two minutes per beam. The PSU sodar integration oc- 


curred over 60 seconds, with each beam sampled every 10 seconds. 
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B. LIDAR THEORY 

The lidar uses the same principal as the sodar to analyze the backscatter from a 
pulse of energv. The lidar utilizes a laser beam to emit pulses of high intensity light that 
will be reflected from gases, aerosols and particulates. A general form of the lidar 


equation is given (Russell et al.,1974): 
DE P, -E AR R) (3.4) 
r p р? , : 


here 

РАК) = received power from range, (R) 
P,=transmitted optical power, 

с = speed of light in air, 

t = duration of laser pulse, 

A = effective telescope aperature, 

B(R) = Optical volume backscattering coefficient and 
1(R) = Optical transmittance of atmosphere. 


This equation is similar in form to the sodar equation (3.1), except that the lidar 
equation is dealing with electromagnetic energy. whereas the sodar equation deals with 
acoustic energy. The received signal from the lidar is a function of the optical properties 
and concentration of the backscattering volume. Since gases and aerosols have different 
effects on the attenuation of light, the received energy can be analyzed to determine what 
causes the scatter. 

Aerosols and gases are distributed by the winds and are trapped in lavers due to the 
stratification of the atmosphere. The detection of these gases and aerosols gives an in- 
dication of the mean wind field and the degree of stratification. The doppler lidar uses 


the same principle of operation to deduce radial velocities as the doppler sodar. 


С. SODAR CROSS SECTIONS 

1. Time-Height Cross Sections 

The wind vector data at the sodar were collected and stored on floppy disk files 

in tabular form, with a header indicating the date and time of the record. The data were 
recorded every 25 meters in the vertical, beginning 50 meters from the surface up to and 
including 730 meters. These files were transferred to the mainframe computer at the 
Naval Postgraduate School on mass-storage for consolidation into separate files of wind 
speed, wind direction, * ertical speed, and o,, 0, and o.. These consolidated files were 
plotted in time-height cross sections. Time for the records was in Universal Time Co- 
ordinates (LTC), but cor verted to decimal days for plotting. The conversion from UTC 
to Pacific Standard Time (PST) is to subtract 8 hours from UTC. Plots were generated 
in the following formats: 24 hours up to 750 meters, four hours up to 350 meters and two 
hours up to 550 meters. The two and four hour plots were centered as close to possible 
to the passage of the sea breeze front. The 24-hour plots started and ended at 0600 local 
time (1400 UTC or 0.5833 decimal day). The raw data from mass-storage were read and 
set to a standaid grid with a resolution of five minutes on the abscissa and 25 meters 
on the ordinate. The density of the grid depended on the time and height resolution of 
the desired field (1.e. 24-hour plots with 750 meter heights were 29 by 288, whereas two 
hour plots with 550 meter heights were 19 by 24). Contours were fitted to the standard 
grid using a rational spline interpolation scheme. Plots were generated and checked for 
bad data, which showed up as unusually strong gradients surrounding a single point. 
Vertical records were manually checked for consistency and errant readings were re- 
moved by comparison to adjacent values. The two hour cross-sections of wind direction, 
wind speed, vertical velocity and c, are included as appendices A, B, C and D. Ihe 
24-hour cross-sections of wind direction and wind speed are included as appendices E 
and E. 
D. DERIVED DATA 

Horizontal divergence and turbulent kinetic energy were derived from the velocity 
ап4 с,. в, and а, fields respectively and plotted in the same format as the velocity fields. 


Horizontal divergence was determined from the continuity equation for 


incompressible flow: 
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The divergence was calculated from the vertical velocity field such that 


Cw җе 
де Уда І. (3.6 
(2 u ) 
This calculation was simple, since the vertical velocity component of the wind with re- 
spect to height was readily available from the sodar data. 


The turbulent kinetic was determined using the values of o,, o, and o, as follows: 
me io? uo? To. (9.7; 
€ = 5 (О, + су ah 1) 


where e equals turbulent kinetic energy. 


E. ANALYSIS OT TIME-HEIGHT CROSS SECTIONS 
1. Wind Direction 
a. Wind Direction Gradients and Onset of Sea Breeze 

The criteria that mark the onset of the sea breeze typically are the change 
in wind direction and increase of wind speed. The wind data from the sodar were ana- 
lyzed to determine the degree of the wind speed change. A sharp change in wind direc- 
tion with tine which indicated sea breeze frontal passage was present in 10 plots during 
LASBEN. Wind direction was contoured with an interval of 30°. Onshore wind was 
Exusidered 270^ + 60° and offshore 90° + 60°. Up valley was 330 + 20°. Down vallev 
was considered 150 + 20°. Frontal passage and the onset of the sea breeze was evident 
when the wind direction showed a strong temporal gradient. The mean temporal gradı- 
ent for sea breezes during LASBEN was on the order of a 90? clockwise shift over 44 
minutes for a strong offshore flow and over 42 minutes for a weak offshore flow. The 
fastest change occurred over 10 minutes during strong offshore flow, while the slowest 
change occurred over 104 minutes during weak offshore flow. 90° was selected as the 
Criterion to define sea breeze onset, when the wind turned through alongshore to 
onshore flow, because the observed temporal wind shift was greatest over this interval. 
The relative strength of the offshore flows was divided into strong and weak flows based 
on the preexisting winds below 400 meters before the onset of the sca breeze. The 
400-meter level was selected since it was the average depth of the established offshore 
flow. Strong flows were considered 3 m's or greater and weak flows less than 3 m s. 
There was no correlation between the sharpness of the front based on the temporal 


gradient and the strength of the previous offshore flow. This agrees with the model 
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result of Pearson et al. (1983) that the intensity of the sea breeze circulation is unaltered 
by the mean flow. Table 2 on page 20 lists the various gradients for sea breeze fronts 
that were detected durıne LASBE® 


Table 2. GRADIENTS OF WIND DIRECTION AT THE SEA BREEZE 
FRONT: p, = mean of strong offshore flows, u, = mean of weak off- 
shore flows 


Temporal Direc- Vert. Directional Offshore Flow (Strong 
tional V(°'min) Vien or Weak) 
ШЕ? 180 37 


ШШ 
e ow 
m ae 
90 42 180 83 па 





A shear laver was evident in the direction plots as a strong gradient in the 
vertical. This shear zone indicates the boundary between the offshore flow and the sea 
breeze. These gradients are shown in Table 2. The mean wind shift in the vertical was 
180* over 38 meters and 83 meters for the strong and weak offshore cases, respectively. 
[t does appear that a strong offshore flow limited the vertical extent of the shear zone 
and that a weak flow allowed a deeper region of directional shear, which agrees with the 
results produced by Simpson and Britter (1980) in a laboratory experiment, as shown in 
Fig. 3 on pape >27, 

In the majority of the frontal passages, the wind direction change preceeded 
a change in wind speed or vertical velocity. In all cases, the direction of the breeze after 
the direction shift was 300° + 20°,which was onshore and up vallev. 

The average onset times for the sea breeze determined at the sodar site by 


the change in wind direction was 1742 UTC. The earliest onset time was at 1644 UTC 
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5, Photographs of a Laboratory Gravity Head Structure in Prevailing 


Flows: The top photo is taken during opposing flow. The bottom is 


taken during accompanying flow. (from Simpson and Britter, 1980) 


and the latest was at 1850 UTC. The time for the increase in wind speed and in vertical 
velocity were also tabulated from the sodar data. 

Onset times were also determined from the surface data. Onset was indi- 
cated by a change in wind direction and speed, temperature and relative humidity. 
Shortly before sea breeze onset, the wind direction was very variable with the trend 


changing from offshore to onshore flow. Wind speeds before onset decreased to near 


zero, then rapidly increased as the sea breeze front passed. Air temperature gradually 
rose due to solar heating in the morning, peaked at frontal passage and then decreased 
as the sea breeze was established. Relative humidity at onset increased quickly. with 
frontal passage and leveled off at the sea breeze continued. A typical surface record with 
sca breeze passage is shown in Fig. 6. 

In most cases, the wind direction shift was consistent with the observed 
sodar wind change. The mean time of onset measured by the surface data was 1813 
UTC. The earliest frontal passage occurred at 1715 UTC. while the latest passage oc- 
curred at 1935 UTC. The variability of wind direction shifts in the surface measurements 
and sodar as well as the absence of data on various davs from both sources explains the 
difference between the surface and sodar mean values. In Table 3, the times of onset 
based on the surface data at the sodar site and the directional shift are shown, as well 


as the times when the speed and vertical velocity increases occurred. 


Table 3. TIMES OF ONSET FOR THE SEA BREEZE: sfcso= Onset based on 
Surface Measurements at the Sodar, sdir= Onset of wind direction shift 
at Sodar, sspd= speed change at Sodar, svrt 2 Time of Maximum Updrafts 
at SOdar XI tines eS 
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Fig. 6. Surface Data Record from 29 September 1987 at the Sodar Site 


b. Height and Duration of the Sca Breeze 

The directional gradient which indicates the boundary between the onshore 
and offshore flow could be followed in the 24-hour contour plots as the sea breeze 
Beenened. ‘ The maximum vertical extent of the sea breeze was determined for the 10 
cases and had a mean elevation of 659 meters, and ranged between 570 and 800 meters. 

The cessation time of the sea breeze was indicated by a counter-clockwise 
turning of the winds, a decrease in the height of the sea breeze and the establishment 
of offshore flow. Duration of the sea breeze was the difference between the onset and 
cessation times determined from the directional data. The mean duration of the sea 
breeze was 10.7 hours. The longest duration was 13.4 hours and the shortest was 8.1 
hours. The mean duration for the sea breeze was also determined using surface data 


measurements at the lidar and sodar sites. Using the surface data, cessation was when 


the wind direction was alongshore and the wind speed had dropped to 2m s. The mean 
durations calculated using surface data at the sodar and lidar were 10.2 and 11.4 hours 
respectively. A summary of maximum height and duration of the sea breeze 15 given in 
Table 4. 


Table 4. MAXIMUM HEIGHT AND DURATION OF THE SEA BREEZE 


Maximum Duration by Duration by Duration by 
Height (meters) Sodar (h) a a Sodar | Surface at Lidar 
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In Fig. 7 on page 25, the increasing height of the sea Breeze ıscelearvsheun as Е 
the maximum height achieved and the decrease in the height with time. 
с. Slope of Sea Breese Fronts 

The slope of the sea breeze front was determined for the first kilometer 
horizontally of the front and for the established sea breeze flow subsequent to the first 
kilometer. The slope for the first kilometer of the front was calculated by measuring the 
difference from the 50 meter level to the highest height of the onshore flow in the first 
kilometer and dividing by the horizontal distance between these points. The slope for 
the rest of the sea breeze was calculated by measuring the difference between the greatest 
height in the first Kilometer and the greatest height of the onshore flow during rest of the 


sea breeze. Horizontal distance was calculated by multiplying the time difference be- 
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24 Hour Directional Plot for 16 September 1987 
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data at the lidar and the time of passage at the sodar. For those cases where a speed 
of advance could hot be calculated, the mean speed of advance was used. The mean 


speed of advance during LASBEX was 2.5 m/s. These speeds are listed in Table 3, 


Table 5. SPEED OF ADVANCE OF THE SEA BREEZE FRONT: Times of 
frontal passage are from the surface data, except in those marked *. which 
are based on sodar data. 


Date Lidar (UTC) Sodar (UTC) A t (h) Speed of Ad- 
vance — 5) 
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In the first kilometer, the mean slope was 1:12 for strong offshore flow and 
1:10 for weak flow. The tvpical siope after the first kilometer was 1:94 for strong off- 
shore flow and 1:160 for weak offshore flow. АП the values except the subsequent slope 
for the weak cases are similiar to ones cited in Simpson et al. (1977), which were 1:4 and 
1:10 for the first Kilometer and 1:100 for the subsequent portion of the sea breeze. 
However, there were no slopes on the order of the one reported by Helmis et al. (1987). 
which had a slope of 1:3 in the first kilometer. There appears to be no correlation be- 
tween the degree of slope in the first kilometer and subsequent to passage. That 1s, a 
strong slope in the first kilometer did not mean that the rest of the sea breeze also had 
a sharp slope, and vice versa. This indicated that the sea breeze was уегу dynamic and 
its velocity and structure changed throughout its duration. The slopes for the sea breeze 
fronts during LASBEX are summarized in Table 6. 


Table 6. SLOPE OF THE SEA BREEZE FRONT: 4, = mean of strong offshore 
flows. u, = mean of weak offshore flows. 


Slope in first kilome- Offshore Flow 

ter meter (Strong or Weak) 
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2. Wind Speed 
a. Synopsis of Wind Speed Data During LASBEX 

A characteristic of the sea breeze that was investigated during this study is 
the difference in time between the change of wind direction and the increase in wind 
speed and whether or not this time difference was affected by the previous offshore flow. 
The general structure of the wind speed field was also compared with the wind direction 
field to see whether features, such as shear zones, indicated in one agreed dimensionally 
with the other, that is if the height and width of the shear zone in one case vas sinillar 
o tlie other. 

The wind speed data were plotted with contour intervals of two meters per 
second for the 24 hour plots and one meter per second for the shorter plots. These plots 
were analyzed for a continuous increase of wind specd with time. The mean time of wind 
Seed increase at the sodar site was 1757 UTC or 17 minutes after the shift in wind di- 
rection. The mean temporal gradient was a 1 m’s change over 55 minutes for strong 
offshore flow and 1 m/s over 21 minutes for a weak offshore flow. In this case, the 
findings do not agree with Pearson's (1983) model that shows no relation between the 


intensity of the mean flow and the intensitv of the sea breeze circulation. The results 


Пот LASBEN point to the shortcommgs that a model has in trying to represent a dv- 
namic flow such as the sea breeze circulation. 

lhere was a shear zone at the boundary between the sea breeze and the 
offshore flow. The vertical gradient of this shear zone was a 1 m.s change over 40 meters 
for the strong case and | m.s over 44 meters for the weak offshore flows. The level of 
the wind speed shear correlated with the level of the wind direction shear. Table 7 lists 


the gradients in speed. 


Table 7. GRADIENTS OF WIND SPEED AT THE SEA BREEZE FRONT: y, 
= mean of strong offshore flow, 4,7 mean of weak offshore flows. 


Vert. Speed (л 00) Offshore Flow (Strong 
c min) or Weak) 
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3. Wind speed afte: frontal passage 

The mean wind speed within three hours after frontal passage during LASBEN 
was 5.3 m s for strong offshore flows and ranged between 4 and 6 m's. The mean for 
weak offshore flow was 7 m/s and was between 6 and 10 m s. The wind speeds of the 
sea breeze reached maximum intensity about the same time as maximum height was 
achieved. The mean maximum speeds during the entire duration of the sea were 8.3 m s 
for stong flows and 7 m's for weak flows. The maximum speeds ranged between 6 and 
10 m's. Offshore flow was clearly present before the onset of the sea breeze. This flow 
was considered the established land breeze and was usually confined below 400 meters. 


The mean speed of the strong offshore flow was 6.3 ms and the mean of the weak 


offshore flows was 2.0 m's. All the weak offshore flows were 2 111 $ or less, while the 
strong offshore flows were between 4 and 8 ms. Table 8 on page 29 gives the various 


maxıma of wind speed. 


Table S. MAXIMUM WIND SPEEDS DURING THE SEA BREEZE AND 
MAXIMUM OFFSHORE FLOW: д, = mean of strong offshore flows. 
f.mu sub w = mean of weak offshore flows. 


Maximum Speed Maximum dur- | Maximum Previ- Offshore Flow 
after onset (m s) ing Sca Breeze ous Offshore (Strong or 
(m s) Flow (m s) Weak) 


10 


а = 
mam; Е 
ЖЕН ШЫ ИШЕ ы (C 1 


Шш. ы 

ЕТІС 2. 
| 
— | ан - 
Први а те | аа 
Sg 7171 1 23 ا‎ ma 





Asthe sea breeze established itself, tne top ofthe onshore breeze increased until 
it assumed a relatively constant elevation. Above the this level, the return flow which 
completes the sea breeze circulation was detectable by the sodar. Since the sodar signal 
had little data past 750 meters, the extent and characteristics of the return flow could 
not be determined. Hlowever, the information from the hdar clearly showed offshore 
flow on 16 September that began around 400 meters and extended vertically up to 1600 
meters. 

4. Vertical Velocity Structure 
a. Vertical Velocity Maxima and Location of the Updraft Region 
The vertical velocity was investigated in order to determine the strength and 
location of the updraft region in the sea breeze front. Vertical velocity cross-sections 
showed good correlation with both the wind direction and speed cross-sections. The 


maximum in upward vertical velocity occurred after the wind direction change and about 


the same time as the increase in wind speed (see Table 3 on page 22). The occurrence 
of the updraft region indicated the front of the longitudinal vortex which is the major 
current within the sea breeze frontal circulation. The mean maximum vertical velocity 
observed after frontal passage was 0.56 m/s, with the range of velocities between 0.2 and 
1.0 m's. This vertical velocity is much smaller than the maximum velocity observed by 
Ogawa et al. (1986) of 3.0 m's; but similar to the values reported by Simpson et al. 
(1977) of 0.7 to 3.0 m/s and the average value that Helmis et al) (1987) observed oie 
m's. A summary of the vertical characteristics of the sea breeze is contain Sii 
Table 9: 


Table $. VERTICAL MOTION CHARACTERISTICS OF THE SEA BREEZE 


Мах. Ver Ез Max. Vert. Ve- Location of Offshore Flow 
tent (m) locity (MS) Max. Velocity (Strong or 
Weak) 
| e om Weak 
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The location of the vertical velocity maximum within the sea breeze 
front gravity head structure is shown in Fig. 8, where the longitudinal vortex turns up- 
ward. 

The location of the updraft region was found to be related to the height of 
the sea breeze (ront, which in turn was related to the strength of the previous oflshore 
flow. The mean location of the maximum vertical velocity was at 207 meters in strong 
offshore flow and 359 meters in weak offshore flows. When offshore winds were strong 


before the sea breeze, the updraft region started at a lower elevation and died out at an 


30 





Fig. 8. Flow regions of the Gravity Head: The stippled region is where the op- 


posing flow is being overrun by the gravity head current (from Simpson 
and Britter.1980) 


average of of 183 meters, and when the winds were weak, the updraft began around 75 
meters and died out around 230 meters. These conditions both agree with the gravity 
head structure in prevailing flow that was observed by Simpson and Britter (1950). In 
a strong opposing current, the gravity head 1s flatter with the flow constrained to a 
shallower region than in onshore or weak offshore flow. The leading edge of the onshore 
flow also tends to have a more pronounced “knee” where it is overriding the counter- 
flow. (see Fig. 5 on page 21) The updraft in the onshore flow is higher, indicating some 
influence from the counter-flow. Fig. S shows how the updraft is affected by the over- 
running and elevated region at the front of the gravity head flow. The “knee” structure 
of the gravity head current was apparent in all the sea breeze fronts investigated during 
LASBEX. 


Sl 


5. Variations of the Vertical Velocity 

The standard deviation of the vertical velocity was used to determine those areas 
of the sea breeze front that were the most turbulent. The standard deviation is higher 
in turbulent regions and near zero in non-turbulent stable lavers. The values for а, were 
directlv generated by the PSU sodar. These values were read and plotted in the same 
manner as the other parameters. Generally, the values of o, decreased with height, 
which indicated a reduction in the amount of turbulence. Vertical velocity fluctuations 
were damped by the boundary at the base of the inversion. Table 10 on page 33 shows 
the change of c, across the front at 80 meters. 

The values of c, did not consistently decrease with time at the lower lerci nie 
< 100 m) after the passage of the sea breeze front. In the majority of the cases, the 
value of e, increased by an average of 0.17 m.s, with other maxima occurring at various 
time in the established sea breeze flow. The changes in o, were between -0.1 and + 0.4 
m.s. Helmis et al. (1987) observed that the values of o, decreased after the passage of 
the sea breeze. whereas. Ogawa et al. (1986). found that the values increased across the 
front. The explanation given by Helmis was that the times of onset and the location are 
the governing factors in the degree of turbulence generated by the sea breeze. In the 
case of LASBEX, the onset of the sea breeze occurred early in the day. The variation 
in vertical velocity was mainly due to mechanical mixing at and within the front. Helmis 
observed sea breezes that passed іп the early afternoon. The sea Dreeze inthisca ea p 
to decrease the amount of thermal activity caused by the heated land. There appeared 
to be no correlation between the strength of the previous offshore flow and the change 
of c, across the front. In the case study for 16 September, the behavior of c, will be 


investigated. 


Е. DERIVED DATA TIME-HEIGHT CROSS SECTIONS 
І. Horizontal Divergence 
Since the sea breeze front is an area of strong convergence, it was logical to 
calculate the intensity of the convergence from the velocity fields. The horizontal di- 
vergence was calculated using the vertical velocity field according to equation (3.6). The 
divergence calculation using the vertical velocity field was straight forward, since the 


sodar produced velocities as a function of z. © 


Table 10. STANDARD DEVIATION OF VERTICAL VELOCITY ACROSS 
THE SES BREEZE FRONT 


o, Before 


After Incrcase or De- | Offshore Flow 
crease (Strong or 
Weak) 
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Strong 
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Lyons and Olsson (1973) calculated and plotted horizontal divergence in time- 
distance coordinates. using the u-component of velocity based on the assumption that 
all variations occurred along a two-dimensional axis (је. Cv/€y = 0), according to the 
continuity equation. The maximum values for convergence that they calculated at the 
ЕШ Гісеге front was on the order of 2 x 10352. The values of convergence calculated 
during LASBEX exceeded J x 102s-! at the head of the front. and at times were as much 
EN 10757. The values calculated from the vertical velocity field of the sodar are 
believed to be a better representation of the actual divergence along sea breeze front. 
because the values of vertical velocity were directly measured. The lake breeze studied 
by Lyons and Olsson had lower computed values of maximum vertical velocity, around 
0.103 то 0.130 га 5, whereas during LASBEX, maximum vertical velocities were almost 
an order of magnitude larger (i.e. 0.8 to 1.0 m's). This is further evidence that Lvons' 
and Olsson’s calculations may have been low and that the assumption of two- 
dimensional flow was perhaps incorrect. Fig. 9 on page 35 shows the typical divergence 
P from Lyons and Olsson. Fig. 10 on page 36 is an example of the horizontal di- 


rgence calculated from the vertical velocity field on 20 September 1987. 


ре 


33 


COMPUTATION 
оао 


A 
« 
~ 
— 
~ 
£ 
M 
x 
9 
ы 
х 


DISTANCE, K^ 





Fig. 9. Divergence values of the Chicago Lake Breeze Study: Contour interval is 


50 х 10:55"), Shaded areas indicate convergence. (from Lyons and Olsson, 
1973) 
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LASBEX HORIZONTAL DIVERGENCE CROSS-SECTION 
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Fig. 10. Divergence values of the LASBEX Study: Convergence at the front oc- 
curs near 20.7071 decimal day. 
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2. Turbulent Kinetic Energy (TKE) 

Since turbulence and mechanical nuxing ıs highest within the head of the sea 
brecze front, the locations of maximum TKE should coincide with those regions of 
strong shear, strong convergence. or overturning flow that creates turbulent eddies. A 
general view of these regions was conducted using the standard deviations from the ve- 
locity fields. The contours for TKE were calculated by squaring the standard deviations 
from the three velocity components obtain their variance for substitution in equation 
(3.7). The values of TKE were then plotted in two hour format. TKE values were found 
to be highest at the base and top of the gravity head, where the longitudinal vortex 1s 
overturning and interacting with both the opposing flow and the vertical shear region 
above the gravity head, corresponding to the area where horizontal flow is directed up- 
wards and into the shear zone within the gravity head. The case study for 16 September 


will show these characteristics. 
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NZ ZCASE STUDWOFR THE SEA BREEZE ON 16 SEPTEMBER 1987 


A. MOTIVATION FOR STUDY 

‘September 16, 1987 proved to be an unusual day since the sodar detected the pas- 
sage of two sea breeze fronts. The sea breeze front that passed in the morning is the 
subject of this case study. The late sea breeze is currently under further investigation. 
Since the 16th was unusual, WPL provided the lidar information gathered during the two 
sea breeze passages. Surface observations and radiosonde observations were not avail- 
able, but were deemed unnecessary in the comparison of the wind field structure of the 


sea breeze that was detected by both the lidar and the sodar. 


B. SYNOPTIC CONDITIONS 

On 16 September, a sea level high pressure center of 1031 mb was located in the 
eastern Nortli Pacific Ocean, northwest of California, along with a weaker high pressure 
Over Arizona. A thermal trough was present in the Central Vallev of California and in- 
tensified during the day, which led to alongshore flow in the afternoon. Winds were 
generaily from the nerthwest along the coast, with a reported wind speed of 5 m’s from 
due west at Monterey in the morning. The wind shifted to the northwesterly in the aft- 
ernoon. Skies were generally clear on the the 16th. The temperature was 56° F at 4 a.m. 
local and rose to a maximum of 76° F later in the day. Fig. 11 on page 38 shows the 


synoptic conditions for this day. 


C. WIND FIELD DATA AND ANALYSIS 
l. Wind Direction and the Onset of the Sea Breeze 

The wind direction cross-section for 16 September was checked for a strong 
change in wind direction with time that indicated the passage of the sea breeze front. 
E се ОГ tle sea breeze occurfed at 1753 UTC (0953 PST). At onset, the wind di- 
теспоп уесгед from 180° to 300°. The temporal wind direction shift in this case was 
90° in 12 minutes. The wind direction cross-section 1s shown in Fig. 13 on page dl. 
Following this directional shift, a perturbation in wind direction occurred directly behind 
the front. The mind diction in this pertubation at 50 meters was 120°, and was con- 
sidered to be the offshore flow being overriden at the base of the sea breeze front. The 
presence of this pertubat:sn may possibly be explained by the “lobe” and “cleft” structure 


that has been observed in the three-dimensional gravity head current. This pertubation 





16 September 


Fig. Il. Synoptic Situation on 16 September: Solid lines are isobars. 


is the area where the opposing current penetrates the gravity current head and causes а 


“cleft” in frontal boundary. The “lobe” is the location of strongest flow at the front of 
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the longitudinal vortex. Ifa “cleft” passed the sodar site, the wind direction would 
change rapidly as the onshore flow of the sea breeze was interrupted by the opposing 
flow and back to the onshore direction as the sea breeze is reestablished. The “lobe” and 
Seem structure 15 shown in Fig. 12. 

A complete directional shift of larger scale occurred about 16.78 decimal day 
(1843 UTC). The wind direction backed with decreasing height, and reached 60° at 50 
meters. This region was at the same location as a maximum in downward vertical ve- 
locity. A possible explanation for this area of directional shift with height is that the 
descending air from the overturning leg of the longitudinal vortex affected the horizontal 
integrity of the flow. The shift in wind direction occurred within the confines of sea 
breeze front which was defined by the directional shear region at its boundaries. The 
structure of the vertical velocity field which complements the wind shift region will be 
discussed later in section 4.3. 

The height of the directional shear zone slowly increased with time. The gradi- 
ent in the this zone was 180° over 83 meters and showed that the layer of marine air 
deepened as the sea breeze became well established. 

2. Wind Speed 

In Fig. 14 on page 42, the increase in wind speed began at 16.75 decimal dav 
(1800 UTC). The wind speed slowly increased from 2 m. s to 3 m.s in about 30 minutes. 
The offshore flow on 16 September was above 8 m s at 450 meters, with speeds greater 
than m s down to 250 meters. When cross-secuons were compared, the change in wind 
speed occurred immediately after the change in wind direction. Helmis (1987) observed 
a wind speed change from 3 m s to 5 m.s in about four minutes across a sea breeze front 
with offshore flow of 5 m s and a change from 2 m, s to 3 m s over 15 minutes with flow 
parallel to the shore. 

Тће speed shear zone dividing the onshore and offshore flow was located be- 
Een 250 and 350 meters. Helmis et al. (1987) found that the top of the inflow laver 
was between 100 and 200 meters; Simpson et al. (1977) observed that the speed shear 
zone was near 600 meters at the head of the sea breeze front and decreased to 300 meters 
after the front had passed; Lyons and Olsson (1973) observed average inflow layer 
heights of 500m, ranging between 100 and 1000 meters; and Ogawa et. al. (1956) ob- 
served that the speed shear layer was between 150 and 320 meters. 

The strong offshore flow (> 4 m/s) was above 250 meters and decreased in in- 


tensity with time as the height of the sea breeze increased. In Fig. 15, the maximum 





Fig. 12. Three Dimensional View of a Gravity Current Head: The “lobe” and 


“cleft structure is shown {топ sonic o 


wind speed after frontal passage was 10 m.s at 500 meters around 17.00 decimal day 
(2400 UTC). The maximum wind speed observed by Simpson et al. (1977) during a sea 
breeze was S m s under 200 meters; by Johnson and O Brien (1973) was 7-5 m Sagem 
elevation of 300 meters; bv Lyons and Olsson (1973) was 6.5 m s at 200 meters; and by 
Ogawa et al. (1986) was 3.5 m s between 100 and 200 meters. The top boundarv of the 
sea breeze flow was not as clearly evident in the wind speed cross-section as it was in the 
wind direction cross-section, but there was a wind speed transition around 550 meters. 
The wind maximum of 8 m.s above this level was associated with the offshore flow. 
3. Vertical Velocity 

In Fig. 16 on page 45, a maximum upward velocity of 0.6 m’s occurred shortly 
after 16.75 decimal day (1800 UTC) and coincided with a “bump” in the wind direction 
contour plot. This “bump” indicated the highest vertical extent of the gravity head 


structure and was present in all the sea breeze fronts that passed the sodar site. 
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Fig. 13. Two-Hour Cross-Section of Wind Direction (16 September 1987): Plot 
Becerra! 7 0°01 C ahrouch 1930 CEC. 


Maximum upward vertical velocities of 3 m s were observed bv Ogawa et al. (1986) and 
of 1.5 m s by Helmis et al. (1987). The top of the vertical velocity maximum, which is 
marked bv the zero vertical velocitv, coincided with the bottom of the shear zone in the 
direction and speed cross-sections. From the direction, speed and vertical velocity 
КО sections; the elevation at the head of the sea breeze front on 16 September was 
between 1S0 and 220 meters. The depth of this same region observed by Simpson et. 


al. (1977) was approximately 600 meters. The updraft region was located between 50 
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Fig. 14. Two-Hour Cross-Section of Wind Speed (16 September 1987) 


and 150 meters. Helmis et al. (1987) observed an updraft region that was from the sur- 
face up te 150) meters: 

A region of downward vertical velocity was present around 16.7917 decimal day 
(1900 UTC). The maximum downward velocity in this region was -0.4 m s. This was 
consistent with descending flow in the rear leg of the longitudinal vortex within a gravity 
head current. Using a mean advance of 2.54 m’s for the sea breeze fronts during 
LASBEX, the horizontal extent of the longitudinal vortex on 16 September was 8.1 km. 


This dimension is similar to the 7 km longitudinal vortex described by Simpson et al. 
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Fig. 15. 24-Hour Cross Section of Wind Speed (16 September 1957) 


(1977). This vortex is an oblong circulation that is located fullv within the the sea breeze 
front. The vortex observed on 16 September was a third as deep as the vortex observed 
by Simpson. This difference between the depths of the vortices was probably due to the 
strength of the offshore flow on 16 September, which was up to $ m's. This strong off 
shore flow flattened the front and restricted the depth of the onshore flow. The de- 
scending leg of the longitudinal vortex observed on 16 September also showed up in the 
wind direction cross section. In Fig. 13. there was a strong temporal gradient of the 


wind direction at 1850 UTC, in the region where the vertical velocity has a gradient of 
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0.4 m.s over 5.6 minutes. This location 1s where vertical motion changed from ascent 
to descent, | 

About 16.7917 decimal day (1900 UTC), the vertical velocity contours indicated 
the area where the sea breeze was reestablished after the penetration of the descending 
leg of the longitudinal vortex. This evidence supports Simpson s theory that the gravity 
head current is a possible description of the actual sea breeze frontal structure. 

4. Standard Deviation of Vertical Velocity 

The location of the strongest fluctuations in the vertical velocity is shown in 
Fig. 17 on page 46. The highest standard deviations of vertical velocity occurred at the 
base of the sea breeze front and around the upward leg of the longitudinal vortex. In 
fact, the whole region of the sea breeze front showed the highest values of standard de- 
viation, indicating strong turbulence and vertical velocity fluctuations in this area. 

The behavior of c, within the boundary laver has been found to follow certain 
relationships (Caughey and Palmer. 1979). Since the c, field was available. this re- 
lationship was investigated. In Caughey and Palmer, c, was normalized using the 
sealing velocity w.. Scaling velocity in the convective boundary layer is determined from 


the heat flux based on the following equation (Ogawa et al.,1986): 

ца (ve, (41) 
where 

g = acceleration due to gravitv (т5 2), 

Ü = mean potential temperature in the boundary layer (°C), 

10”, = sensible heat flv“ at surface (ms? Cand 

z; 2 depth of the mixed laver (72). 


The value of w., could not be calculated since heat flux data were not available. How- 
ever, over a period 50 minutes, the change in heat flux was considered to be negligible 
and, therefore, the value of w. was also taken to be constant. The value of z, was found 
as the height on the o, cross-section which was indicative of the base of the inversion 


(1.е. where the values are decreasing quickly towards 0). 
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Fig. 16. Two-Hour Cross-Section of Vertical Velocity (16 September 1987) 


The values of o, were plotted versus z/z, to compare with Ogawa et al. (1956). 

he emphasis of this comparison was on the distribution of o,. In Fig. 18 on page 47 
and Fig. 19 on page 4S, the maximum value in o, occurred at z/z,= 0.25 for the pre- 
frontal case and at z/z,= 0.3 for the postfrontal. The heat flux was allowed to vary, so 
that the maximum in Caughey and Palmer’s (1979) curve matched the variance of the 
acoustic sounder values of o, are at a maximum in comparison with z/z, . The plot by 


Ogawa et al. (1986) is shown in Fig. 20 on page 49. 
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Но. 17. Two-Hour Cross-Section - Standard Deviation of w (16 September 1957) 


Ogawa observed that the maximum value of the normalized o, occurred at z/z, 
= 0.3. The values of c, decrease with height as Fig. 20, which contains the free con- 
vection case as described by Caughey and Palmer (1979). Contrary to Ogawa, the values 
of c, in the LASBEX plots did not fall off as dramatically with height and at times ac- 
tuallv increased with height. This was because of the strong offshore flow above 400 
meters on 16 September caused more turbulence above the inversion. The wind profile 
that Ogawa observed decreased in intensity with height above the inversion, which gave 


lower values olo" 
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Fig. 18. . Standard Deviation of w on 16 September: co, versus z/z, for the Pre- 


frontal Case. Solid curve = Caughey and Palmer (1979) 
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Fig. 19. Standard Deviation of w on 16 September: o, versus 2/2, for the 


Postfrontal Case. Solid curve = Caughey and Palmer (1979) 
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Fig. 20. Standard Deviation of w normalized: с, normalized versus 2/2. Solid 


curve = Caughev and Palmer (1979). (from Ogawa, 1986) . 
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Horizontal Divereence 

The divergence cross-section for the 16 September is shown in Fig. 21 on page 
51. When compared with the wind direction cross-sections, the horizontal divergence 
cross-section showed an area of convergence at the head of sea breeze front. Fig. 21 
shows that the maximum values of convergence occurred at the head of and within the 
sea breeze front. Values of convergence in these regions were up to Š x 102s-!, The first 
region is located at the top of the maximum in Vertical velocity and the second 1s located 
where the the longitudinal vortex is turning back into the foward flow. What is partic- 
ularly notable ıs that this cross-section showed that the convergence followed the wind 
shift boundary of the sea breeze front outlining the boundary of the gravity head struc- 


ture, similar to the divergence fields of Lvons and Olsson (1973). 
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Fig. 21. Two-Hour Cross-Section of Horizontal Divergence (16 September 1957) 


6. Turbulent Kinetic Energy 
The cross-section ofthe TKE for 16 September showed that the highest values 
were located ın the sea breeze front where mechanical mixing and vertical motion gen- 
erated the greatest turbulence. The actual maxıma in turbulent kinetic energy were lo- 
cated in the region where the overturning air turned back into the forward flow and at 
the strong shear zone near the base of the front. The TKE cross-section for this day is 


shown in Figs 22 om pace 52: 
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Lim». Two-Hour Cross-Section of Turbulent kinetic Energy (16 September 
1987) 


D. LIDAR DATA 
l. Color Coding 
Color coding was used on the WPL lidar plots to show the intensity and direc- 
tion of the radial velocity field. "Cool" colors were used for velocities towards the lidar, 
Whereas “warm” colors were used for velocities away from the site. The colors in each 
case were graduated with the “cool” velocities graduated from green to blue to purple, 


with green indicating slower speeds and purple indicating the highest speeds. The 
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graduation for the “warm” colors was from tan to yellow to red, with tan being the 
slowest speeds and red the highest speeds. White was used to indicated zero radial ve- 
locities. 

2. RHI Plots 

The lidar information available for 16 September 1987 closely agreed with the 

wind field produced by the sodar. The RHI on this day was primarily oriented on a line 
of 90°. The separation between the low level sea breeze and the intense offshore flow 
was Clearly indicated. The depth of this offshore flow started at 500 meters and had a 
depth of 1200 meters. The flow above 2 km was onshore. About 1722 UTC, there was 
an indication that the sea breeze front was moving away from the lidar site. This is 
shown in Fig. 23 as an increase in the speed of the flow moving away from the lidar to 
the east (right side of picture) at the lowest levels. 
By 1735 UTC, the sea breeze had an estimated depth of about 200 meters, based on the 
height of the radial velocity field in the lowest layer. This agreed with the depth of the 
onshore flow determined bv the sodar, which was between 200 апа 250 mete ые 
actual location of the sea breeze front was not determined due to the spatial resolution 
of the output information. 

There was one RHI view of the sea breeze oriented down the Salinas Valley of 
140&degree; it is shown in Fig. 24. The onshore flow was easier to detect at this ori- 
entation because it 
showed up as a deeper and stronger flow both towards and awav from the lidar. The 
speed of the offshore flow from the Salinas Valley was weaker, with a radial velocity of 
dms. The strong flow above the sea breeze (> 6 ms) was believed to be due tome 
geography of the area and is explained in the next section 

3, PPI Plots 

The PPI plots give a better indication of the offshore flow region and the 
progress of the sea breeze front. The 5° clevation angle plots show a velocity maximum 
at 5320 meters of 8 m/s. This is shown in Fig. 25 as a velocity maximum towards the 
lidar from the northeast and the corresponding maximum away from the lidar towards 
the southwest. The elevation and orientation of this radial velocity maximum corre- 
sponded very well with the sodar wind direction and wind speed cross-sections. A region 
of wind speed of 7 m's with a direction between 60° and 90° was present above the sodar 
site at the same time as the radial velocity maximum was above the lidar site. (Refer to 


Fig. 13 on page 41 for the wind direction and Fig. 14 on page 42 for the wind specd.) 
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Gece iit Lidar Display at 1722 UTC, 16 September 1937: Top picture 1s the 
REI display. The lidar is at the center. Due east is to the right side of the 


picture. Distances and elevations are in kilometers. Velocities are in m s. 
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Fig. 24. RHI Lidar Display at 1805 UTC, 16 September 1987: Top picture is 
RHI display. Down vallev is to the right of the picture. Distances and el- 


evations are in kin. Velocities are in m s. 
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Fig. 25. PPI Lidar Display at 1716 UTC, 16 September 1987: Elevation Angle 
was 5°. True north is at the top of the picture. Rings are spaced 4 kilo- 


meters apart. Velocities are in m.s. 
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According to the lidar, the offshore flow was from the high terrain to the 
northeast and from the gap in the coastal ranges where the San Juan Valley ts located. 
The strength of this offshore flow was attributed to the fact that the San Juan Vallev is 
an outlet of airflow for the higher Santa Clara Valley which induced a strong offshore 
current. 

The PPI piots at 1° elevation showed the progression of the sea breeze front. 
The convergence at the head of the front was detected as it moved in a general southeast 
direction away from the lidar site. The lidar showed that the sea breeze front was in- 
muencedsox the surrounding terrain. At 1719, the sea breeze front was estimated to be 
located at or around the lidar site, since some weak convergence was approximately 1.5 
km south. As shown in Fig. 26. the sea breeze flow advanced northeast up to the high 
terrain approximately 7 km due east of the lidar and north to the mouth of Elkhorn 
Slough. At this time, the sea breeze front had not advanced at all into the Salinas Vallev. 

The high terrain to the east effectively blocked the further advance of the sea 
breeze front at lower levels. The calculated height of the terram at 7 km was about 120 
meters. By 1743, the front had moved about 4.8 km south and southeast from the lidar 
and had began to move up Elkhorn Slough, north of the lidar. This is shown in 


Fig. 


а 


27 on page 58. 

The estimated time of sea breeze frontal passage at the sodar site based upon 
the lidar data was around 1800 UTC and compared favorablv with the time determined 
from the sodar wind field data . which was 1745 UTC. Using the lidar positions of the 
convergence zone at the sea breeze front. the calculated speed of advance was 2.09 m s. 
This compared well with average speed of 2.54 m’s calculated from surface data. 
Fig. 25 on page 39 indicated that the sea breeze had completely filled the mouth of the 


Salinas Valley and all low terrain to the east of the lidar by 1827 UTC. 
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PPI Lidar Display for 1719 UTC. 16 September: Elevation Angle was 
1°. True north is at the top of the picture. Rings are spaced 4 kilometers 


apart. Velocities are m п s. 
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іс. 27. PPI Lidar Display for 1747 UTC, 16 September 1987: Elevation Angle 


was 1°. True north is at the top of the picture. Rings are spaced 4 kilo- 


meters apart. Velocities are in m s. 
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Fig. 28. PPI Lidar Display for 1827 UTC. 16 September 1987: Elevation Angle 


was 1°. True north is at the top of the picture. Rings are spaced 4 kilo- 


meters apart. Velocities are in ms. 
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V. SUMMARY AND CONCLUSIONS 


A. MAPPING THE SEA BREEZE CIRCULATION 
l. Accuracy of the Sodar and Lidar 

Various features that have been produced in laboratory experiments such as the 
flow regions of a gravity head current were observed during LASBEX, most notably the 
structure 1n the frontal region where turbulence and mechanical mixing produce large 
Variations in velocity. This area is caused by the overturning of the longitudinal vortex 
that makes up the major flow within the frontal region and its interaction with the am- 
bient flows surrounding it. 

The unique combination of the lidar and the sodar proved to be an effective 
method to detail the structure of the sea breeze circulation. The instruments both per- 
formed well in a marine environment. The sodar deteeted and mapped the passage of 
the sea breeze front, whereas the hdar showed the separation between the onshore and 
offshore flow. as well as the progress of the convergence zone associated with the sea 
breeze front. 

2. General Characteristics of the Sea Breeze Circulation 

The general characteristics of the sea breezes that were observed during 
LASBEX on the Monterey Bay Salinas Valley in September 1987 were similar in manv 
Ways to those sea breezes seen in other studies, 

The wind direction change over time was found to be an accurate indicator of 
sea breeze onset." The times of onset were consistent, occurring around 1742 UTC at the 
sodar station. which was 6.5 km inland. Differences in the relative slope, speed and 
general shape of the sea breeze flow was observed for each front as it passed the sodar. 
The directional shear in the vertical meausred by the sodar gave a good indication of the 
depth of the onshore flow. The lidar presented a similar depth in the RHI scan. 

The observed sea breezes advanced onto land with a mean speed of 2.5] m's. 
and lasted between 8 and 13 hours. The maximum height that the onshore flow achieved 
averaged about 659 met ‘rs vertically. The slope of the onshore flow measured using the 
height of the directional shear zone gave values that ranged between 1:5 and 1:17 in the 
first kilometer horizontaiiy of the flow and between 1:56 and 1:220 in the established 
onshore flow. These values for the slope of the onshore flow compare favorably with 


ese cited by Simpson et al. (1977). 
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The increase of wind speed and the maximum updraft region both occurred after 
the change in wind direction. This agreed well with the location of the longitudinal 
vortex within the sea breeze front. The maximum wind speeds of the sea breeze recorded 
at the sodar were from 6 to 10 ms. Vertical velocities observed during the sea breeze 
in the mouth of the Salinas Valley ranged between 0.2 and 1.0 m's and agreed with 
typical values recorded in previous studies, such as Ogawa et al. (1986). The correlation 
between depth of onshore flow and the vertical extent of the updraft region was ob- 
served. A larger area of upward flow was associated with weak offshore flow and a 
deeper front. The location of downdraft regions corresponded to the descending leg of 
the longitudinal vortex and gave a rough measure of the length of this current. On 16 
September, the length of such a flow was found to be around 7 km long, and compares 
mellowith tne vortex observed by Simpson et al. (1977). 

The horizontal divergence calculations made using the vertical velocity field 
pusuuced values of convergence that at times were up to 12x10s7, These were found 
to be twice as large as the values calculated bv Lyons and Olsson (1973) in a Chicago 
ІП бгееле study, The convergence values calculated from the LASBEX vertical veloc- 
itv field are believed to be the most accurate for the sea breeze front to date. 

3. Tracking of the Sea Breeze Front 
The lidar was able to detect and track the convergence zone associated with the leading 
edge of the sea breeze front. The derived mean speed of advance agreed favorabh with 
speeds derived through surface measurements at both the hdar and sodar, as well with 
speeds seen in other studies. "The depth of the offshore flow detected by the lidar was 
over twice as large as the onshore flow, which is the typical dimension seen in other 
studies. However. the speeds in this offshore flow were either the same of larger than the 
onsnore flow. This observation is contrary to those studies which saw that the speed 
ol the oilshore flow was roughly half as fast as the onshore flow." In PPI, the lidar 
showed how terrain affected the movement of the sea breeze front, such as halting the 
advance of the sea breeze front and the movement of the onshore flow into Elkhorn ` 
Slough. 


B. RECOMMENDATIONS 
Since many of the features of the sea breeze flow have been observed by the sodar, 
further studies should concentrate using this versatile instrument to complete the dv- 


namic structure of the sca breeze front, with concentration on the variations of the 


61 


velocity components. These variations would help understand the dimensions and 
structure of the turbulence within the sea breeze front. 

As far as improvements in the mapping of the sea breeze structure, there was an- 
other sodar which operated concurrently with the PSU sodar, but whose data was not 
available at the time of this study. The combination of this study and the data from the 
other sodar would give a better view of the sea breeze structure and refine the actual 


spatial dimensions of its associated flows. 


APPENDIX A. TWO-HOUR WIND DIRECTION CROSS SECTIONS 
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Fig. 29. Two-Hour Wind Direction Cross Section - 17 September 1987 
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ig. 30. Two-Hour Wind Direction Cross Section - 18 September 1987 
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Two-Hour Wind Direction Cross Section - 19 September 1987 
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Two-Hour Wind Direction Cross Section - 20 September 1987 
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Two-Hour Wind Direction Cross Section - 22 September 1987 
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Two-Hour Wind Direction Cross Section - 29 September 1987 
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Two-Hour Wind Speed Cross Section - 17 September 1987 


Fig. 38. 
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Two-Hour Wind Speed Cross Section - 18 September 1987 
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Two-Hour Wind Speed Cross Section - 19 September 1987 
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Two-Hour Wind Speed Cross Section - 20 September 1987 


Fig. 41. 
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Two-Hour Wind Speed Cross Section - 21 September 1987 
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Two-Hour Wind Speed Cross Section - 23 September 1987 


Fig. 44. 
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Two-Hour Wind Speed Cross Section - 27 September 1987 
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Two-Hour Wind Speed Cross Section - 29 September 1987 


Fig. 46. 
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Fig. 47.  Two-Hour Vertical Velocity Cross Section - 17 September 1937 
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Two-Hour Vertical Velocity Cross Section - 21 September 1987 
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Fig. 53.  Two-Hour Vertical Velocity Cross Section - 23 September 1987 
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Two-Hour Standard Deviation Cross Section - 18 September 1987 
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Two-Hour Standard Deviation Cross Section - 19 September 1987 
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Two-Hour Standard Deviation Cross Section - 22 September 1987 
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Two-Hour Standard Deviation Cross Section - 29 September 1987 
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A. NOTES 
The cross-section for 23 September 1987 was too noisy between 23.958 and 24.208 
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Fig. 65. 24-Hour Wind Direction Cross Section - 17 September 1987 
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24-Hour Wind Direction Cross Section - 18 September 1987 
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24-Hour Wind Direction Cross Section - 19 September 1987 
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Fig. 68. 24-Hour Wind Direction Cross Section - 20 September 1987 
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Но. 69. 24-Hour Wind Direction Cross Section - 21 September 1987 
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Fig. 70. 24-Hour Wind Direction Cross Section - 22 September 1987 
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Fig. 72. 24-Hour Wind Direction Cross Section - 27 September 1987 
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Fig. 73.  24-Hour Wind Direction Cross Section - 29 September 1987 
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24-Hour Wind Speed Cross Section - 17 September 1987 
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Fig. 75. 24-Hour Wind Speed Cross Section - 18 September 1987 
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Fig. 76.  24-Hour Wind Speed Cross Section - 19 September 1987 


116 


03$/SU3L3N 2 :IVAU3LNI UNOLNOD 


(Ava 1VWIO3Q)3WIL 
6051 99912 сее іе 803123 ево іе 858702 сев'0г 80202 


' 
' 
.. 
' 
' 
' 
' 
' 
о 
а 


¿861 HAEWALdAS 
МОІ1935 55099 A33d4S (QNIM X38SV'1 


68502 


0:059 0055 05» 0056 0052 005: 0706 
(SY3L3W) LHOI3H 


0'052 





24-Hour Wind Speed Cross Section - 20 September 1987 


piu 


Fig. 


117 


O3S/SUu313N 2 -IVAU3LNI UNOLNOD 


(AVG 1VWIO3O)3WIL 
ees 32 esp zz 666722 802726 68022 868712 668% 807712 





% 
4 
>10 
0 

! 
% 
, 

` 
` 


486. H38I/N31d3S 
NOILO3S SSOUuO Q33d$ QNIM X38SV1 


68512 





0'059 0'0SS 0705» 0056 0053 07061 0705 
(SH313W) LHOI3H 


0'054 


118 


24-Hour Wind Speed Cross Section - 21 September 1987 
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